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Abstract 

Using data taken during the year 2000 with the NA48 detector at the CERN SPS, 

a search for the CP violating decay Kg Stt^ has been performed. From a fit to 
the hfetime distribution of about 4.9 milhon reconstructed / ?,tt^ decays, 
the CP violating amplitude 77000 = ^(-^^s ^t^^) / A{Ki^ Stt'^) has been found 
to be Re(r/ooo) = -0.002 ± 0.011 ± 0.015 and Im(r/ooo) = -0.003 ± 0.013 ± 0.017. 
This corresponds to an upper limit on the branching fraction of Br(K^s Svr^) < 
7.4 X 10~'^ at 90% confidence level. The result is used to improve knowledge of Re(e) 
and the CPT violating quantity Im(5) via the Bell-Steinberger relation. 
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1 Introduction 



The violation of CP symmetry was discovered in 1964 in the decay of the long-hved 
meson to two charged pions [1]. Since then, other CP violating Ki^ decay modes 
— in particular direct CP violation — and CP violation in decays have been 
observed. The short- hved neutral kaon Xg, too, should manifest CP violating decay 
amplitudes. However, due to the large width of the meson, the branching ratios of 
CP violating decays are 6 orders of magnitude smaller than the branching ratios 
of the main CP violating Ki^ decays. One unambiguous signature of CP violation in 
Kg decays would be the observation of the decay — > Stt". Within the Standard 
Model its branching ratio is predicted to be 1.9 x 10~^. 

Since the kaon is spinless, the Stt" final state has a well-defined CP eigenvalue of 
— 1, when neglecting direct CP violation, and the decay K^, — > 37r° is CP forbidden. 
The CP violating parameter r^ooo is defined as the amplitude ratio 

_ A{K, ^ 37r°) 



When assuming CPT conservation and neglecting isospin 7 = 3 and non-symmetric 
1=1 final states, r/ooo = ^ + i Im(Ai)/Re(Ai) with the parameter e of indirect CP 
violation and the 1=1 amplitude Ai. The imaginary part of r/ooo is in principle 
sensitive to direct CP violation [2]. 

Additional interest in — > 37r'^ decays arises from the search for CPT violation. 
The Bell-Steinberger relation [3] links possible CPT violation in the K^K^ mixing 
matrix with CP violating amplitudes in and Kc, decays via the conservation of 
probability. At present, the hmit on CPT violation is limited by poor knowledge of 

^000- 

In a recent investigation, the CPLEAR experiment found Re(?7ooo) = 0.18 ± 0.15 
and Im(?7ooo) = 0.15 ± 0.20 in the decay of flavour-tagged and mesons [4], 
which corresponds to an upper bound of Br(irg — > 37r°) < 1.9 x 10~^ at 90% 
confidence level. In addition, the SND collaboration has set a limit of Br(irg — > 
37r°) < 1.4 X 10-5 [5]. 
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Fig. 1. Beam- line of the NA48 experiment during the 2000 data taking. 



In this Letter we report an improved measurement of the r^ooo parameter using data 
collected from a short neutral beam with the NA48 detector at CERN. Sensitivity 
to 7^000 comes from Kq/K^ — >■ Stt" interference at small decay times near the target. 
A beam of pure K^^ decays, taken under the same experimental conditions, was 
used for normalization. The experimental set-up is described in the next section, 
while section 3 deals with event selection and reconstruction. The fit to extract the 
parameter t^qoo is described in Section 4. The last section discusses the implications 
of the result with respect to the 37r° branching fraction and the limit on CPT 

violation. 



2 Experimental Setup 

The NA48 experiment was designed for the measurement of direct CP violation in 
neutral kaon decays. Its main feature is two simultaneous, almost collinear beams 
of neutral kaons derived from proton beams from the CERN SPS delivered to two 
fixed targets [6]. The kaon beams have a common decay region and decays from 
both beams are recorded with the same detector (see Fig. 1). Both targets are made 
from beryllium and have a length of 400 mm and a diameter of 2 mm. The far 
target is located 126 m before the beginning of the decay region while the near 
target is only 6 m up-stream of the fiducial region and displaced by 7.2 cm in the 
vertical direction from the axis of the far-target beam. The two beam axes have an 
angle of 0.6 mrad with respect to each other and cross at the longitudinal position 
of the electromagnetic calorimeter, 120 m down-stream of the near target. In both 
beams, charged particles are deflected by sweeping magnets. Particle decays from 
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the far target are almost exclusively K^^ decays, while decays originating from the 
near target are mainly Kg (and neutral hyperon) decays with, however, a small 
component of K^^ decays. 

The analysis reported here is based on data from a 40-day run period in 2000, with 
only the near target in operation. The neutral beam was produced by 400 GeV/c 
protons at a production angle of 3.0 mrad, with a beam intensity of about lO^*' 
protons during a 3.2 s long SPS spill. This was about a factor of 300 higher than 
the typical intensity of the near-target beam during the direct CP violation mea- 
surement. This high-intensity run period was preceded by a 30-day run with only a 
far-target beam under the same beam conditions, except for the proton beam energy 
being 450 GeV and the production angle being 2.4 mrad. The data from this first 
run period were used for normalization purposes. 

In the usual configuration, the main NA48 detector elements are a magnetic spec- 
trometer for reconstruction of charged particles, followed by a hodoscope for charged 
particles and a liquid krypton electromagnetic calorimeter (LKr). During both run 
periods in 2000, the spectrometer drift chambers were absent and its vessel evacu- 
ated, leaving no material between the final collimator and the hodoscope for charged 
particles directly in front of the calorimeter. In the far-target run period, the spec- 
trometer magnet was powered with an integrated field of 0.88 Tm in order to have 
the same running conditions for systematic studies of the direct CP violation mea- 
surement. In the subsequent near-target run period the magnet was off. 

The liquid krypton calorimeter measures the energies, positions, and times of electro- 
magnetic showers initiated by photons and electrons [7]. It has a length of 127 cm, 
corresponding to 27 radiation lengths, and consists of 13212 cells in a projective 
tower geometry which points to the middle of the decay volume. The active volume 
has an octagonal shaped cross section of about 5.5 m^. Each cell has a 2 x 2 cm^ 
cross section and is formed by copper-beryllium ribbons which are extended longi- 
tudinally in a ±48 mrad accordion structure. The cells are contained in a cryostat 
filled with about 10 m^ of liquid krypton at a temperature of 120 K. The initial 
ionization signal induced on the electrodes is amplified, shaped, and digitized by 
40 MHz FADCs. The energy resolution of the calorimeter was 

a(E) 0.090 0.032 



with E in GeV. The spatial and time resolutions were better than 1.3 mm and 300 ps, 
respectively, for a photon with energy above 20 GeV. The read-out system was 
calibrated by a charge pulse every burst during data taking. The relative calibration 
of the individual cells was determined by using K^s decays during the 1998 run period 
and checked to be similar in 2000 by n^/i] 77 decays produced in thin plastic 
targets in a special run with a tt^ beam. The final calibration of the overall energy 
scale was performed by fitting the effective edge of the coUimator as reconstructed 
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in the data to the one in the Monte Carlo simulation. 

The hadron calorimeter follows the LKr calorimeter and was used as a veto-counter. 
It is composed of forty-eight 24 mm thick steel plates interleaved with scintilla- 
tor plates and measures energies and horizontal and vertical positions of hadronic 
showers. A more complete description of the NA48 detector can be found in Ref . [8] . 

The trigger decision for 3n^ events was based on projections of the deposited energy 
in the liquid krypton calorimeter [9]. In both run periods the trigger required a 
total deposited energy of at least 50 GeV. In addition, the radius of the energy 
centre-of-gravity had to be less than 15 cm from the detector axis and the proper 
kaon lifetime, measured from the final collimator, had to be less than 9 Kg lifetimes 
with both radius and proper lifetime computed online from the moments of the 
energy depositions in the calorimeter. The trigger efficiency was determined using 
a 37r° data sample triggered by a scintillating fiber hodoscope located inside the 
calorimeter. It was measured to be 99.8% in both the near- and the far-target run 
periods, and showed no dependence on energy or decay vertex position within the 
decay volume used in the analysis. 



3 Event Selection and Reconstruction 

To identify 871° — * 67 events and to determine their kinematics, the measured 

energies positions Xi and i/j, and times of the photon showers in the liquid krypton 
electromagnetic calorimeter are used. Prom the energies deposited in the LKr cells, 
clusters were formed, which had to fulfill the following selection criteria. The cluster 
energies were required to be above 3 GeV and below 100 GcV, within the range of 
the linear energy response of the calorimeter. To avoid energy losses, each cluster 
had to be more than 5 cm from the edge of the beam pipe and from the outer edge of 
the sensitive area. In addition, the distance to the closest dead cell of the calorimeter 
was required to be larger than 2 cm. 

On all combinations of 6 clusters which passed these requirements, the following 
further selection criteria were applied. To avoid difficulties with the sharing of energy 
among close clusters, the minimum distance between two clusters had to be at least 
10 cm. All six clusters were required to lie within 2 ns of the average cluster time. 
The sum of cluster energies had to lie above 60 GeV and below 185 GeV. The radial 
position r^og of the energy centre-of-gravity at the calorimeter had to be less than 
7 cm in the near-target and less than 4 cm in the far-target run, which had stronger 
beam coUimation. 

To avoid background from accidental pile-up, events with more than one such com- 
bination of clusters and those with additional clusters with an energy above 1.5 GeV 
within 3 ns of the event time were rejected. 
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Fig. 2. Longitudinal vertex position of reconstructed — > 3tt^ events. In the figure, the 
far-target data have been normaUzed to the near-target data in the fiducial decay region. 

From the surviving candidates, the longitudinal vertex along the z direction was 
reconstructed from all pairings by assuming the decay of a with the nominal 
mass rriK- 
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i=l j>i 



(3) 



with the distance Zi^j^^ between the near target and the calorimeter. The resolution 
of the reconstructed longitudinal vertex position is about 60 cm, corresponding to 
~ 0.1 Kg lifetimes for typical kaon energies. 

Using the longitudinal vertex position ;Zvcrtcx; the invariant two-photon masses mi, 
7712, and 777,3 were computed for all 15 possible photon pairing combinations, and a 
X^-like variable was constructed as 

2 /i(mi + m2 + ms) - m,,i\" /i(mi - ™±™)y / i(,„2 _ majV ,„> 
= ^ a, ^ — 



where the resolutions o"j, parameterized as a function of the smallest photon energy, 
were determined from data. The combination with the lowest value of Xs-kO was 
chosen. In addition, a value of xln^ l^ss than 90 was required for this combination. 

To reject possible background from hadrons, we required the total energy deposited 
in the hadron calorimeter within 15 ns of the event time to be less than 3 GeV, 
which removed 0.15% of the signal events. 

The -Zvortox and energy distributions of the selected events are shown in Figs. 2 and 
3. For the fit of the parameter 77000, only events with a longitudinal vertex position 
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Fig. 3. Energy spectra of reconstructed 3tt^ events. The far-target data have been 

normahzed to the near-target data. 
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Fig. 4. Energy versus hfetime, measured from the target, of accepted K Svr events 
from the near-target run. The contour encloses the accepted region for the ryooo fit- 

^vertex > 8 m were considered to avoid detector resolution effects for vertex positions 
near the final collimator at 2; = 6 m. The down-stream vertex region was limited by 
^^vcrtcx < 55 m and a maximum lifetime of Sr^g- In addition, events with lifetimes 
tcou/^Ks > 0-8 + 0-06 X Ex/GeV, measured from the final collimator, were rejected 
to avoid a region at low energies and high lifetimes where the trigger was partly 
inefficient. The accepted region in proper time and energy is shown in Fig. 4. 

In total, about 4.9 x 10® K^^ s ~^ Svr^ events were reconstructed from the data of the 
near-target run and about 109 x 10® Kj^ Stt® events from the far-target run data. 



10 



4 Data Analysis 



4.1 Method of the Measurement 

At the targets, and mesons are produced by strong interactions in equal 
amounts. The K — > 37r° intensity as a function of proper time t, measured from the 
target, is then given by 



with the total Kj^ and Kg widths Fl and Fs and the Kj^Kg mass difference Am. 
The dilution D{p) — {Nko — Nj^)/{Nko + N-^) describes the momentum dependent 
production asymmetry between and at the target. For the 77000 measurement 
we analyzed data from the near-target run period, using the pure K^^ 37r° from 
the immediately preceding far-target run period to correct for trigger, acceptance, 
and reconstruction efficiencies. The difference of the kaon beam momentum spectra 
between the two periods was taken into account by performing the analysis in 5 GeV 
wide bins of energy covering a range from 60 to 185 GeV. The two set-ups had small 
differences in geometry, incident kaon beam angles, and collimation. Based on Monte 
Carlo studies with samples of about 90 million reconstructed events for each beam, 
the ratio 



with the numbers A^°'=^" and A^'" of reconstructed events and the acceptances e"""' and 
gfar fQj. j^gg^j. £g^j. targets was determined. The geometrical correction e"'''"(t)/e*'''(t) 
between the two beams nearly equals to 1 and is shown in Fig. 5 for three different 
energy intervals. Neglecting this correction would shift the 77000 result by ARe(77ooo) = 
-0.03 and AIm(77ooo) = 0.03. 

The intensity /a^o (t) as function of proper time might be altered if there is significant 
Ks regeneration by K^^ mesons hitting the final collimator. This possibility has been 
considered and found to be negligible for the measurement presented here. 

In the near-target data decays K — > tt'^tt^tt^ with one pion undergoing a Dalitz decay 
TT^ e+e~7 had to be taken into account. Since one shower could have missed the 
detector or could have had an energy below the detection limit of 1.5 GeV, these de- 
cays could be mis-identified as K ^ Sir^ 67 events. In many such cases, due to the 
energy loss, the decay vertex position z^^,^^^ and the kaon lifetime were reconstructed 
further down-stream (see Eq. 3). For the near-target run the mis-identification rate 



/3^o(i) oc e 




(5) 




(6) 
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Fig. 5. Ratio g^'^^'^/gf^'^ of the acceptances for 3tt^ events from the two kaon beams in three 
kaon energy ranges, as determined by Monte Carlo simulation. 



with respect to the acceptance of good 3n^ events was e(27r°7r£))/e(37r°) ^ 40%, 
while it was neghgible in the far-target run due to the presence of the magnetic 
field. Dalitz decays were taken into account in the Monte Carlo generation of the 
near-target K^^ Sir^ events. 



4-2 Fit to the data 



For all 37r° candidates, the proper lifetime of the kaon, measured from the po- 

sition of the near target at 2; = 0, is computed as t = -Zvertox/ (7/3c) ~ ^^vertex ^k/ (Ekc), 
where mx is the kaon mass taken from the PDG and is the sum of the cluster 
energies. 
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Fig. 6. The ratio of near-target over far-target Stt^ data, corrected for acceptance differ- 
ences, for three different energy intervals. The points with error bars are the data. The 
continuous curve, which practically coincides with 1, shows the result of the fit for 77000- 



A global least-squares fit to the time evolution /3^o(i(:) (Eq. 5) is performed on the 
corrected time distributions fsnoit) (Eq. 6) for each kaon energy interval. Free pa- 
rameters in the fit are Re(?7ooo)5 Ini(^ooo), and the normalization constants of each 
energy interval. The K^K^ dilution is assumed to be linearly dependent on energy, 
as expected from the parton model. It is obtained using data from the NA31 experi- 
ment [10], corrected for the different momentum spectrum and production angle [11] 
and parameterized as a linear function of the energy. 



The result of the fit gives Re(r/ooo) = -0.002 ±0.011 and Im(r7ooo) = -0.003 ±0.013 
with a correlation coefficient of 0.77 between real and imaginary parts (Fig. 6). The 
of the fit is 689 with 660 degrees of freedom. 
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4-3 Systematic uncertainties 



Wc have investigated contributions to the systematic uncertainties from detector 
acceptance, from calorimeter energy scale and non-linearities, from possible back- 
grounds and accidental activity, and from K^K^ dilution. 

As described above, the acceptance determination relied on data from the far-target 
run, with the Monte Carlo simulation used to calculate residual geometrical differ- 
ences. Many checks of the acceptance calculation were performed. Differences in the 
reconstruction efEciency between the two runs were investigated by varying the Xg^o 
cut, which led to uncertainties of ARe(77ooo) = ±0.009 and AIm(r7ooo) = ±0.013. 
Remaining resolution effects near the collimator were studied by varying the up- 
stream vertex cut and gave an uncertainty of ARe(ryooo) = ±0.010, while the imagi- 
nary part — playing a role only at larger lifetimes — was unaffected. By comparing 
two Monte Carlo samples with different target-coUimator geometries, we estimated 
ARe(?7ooo) = ±0.005 and AIm(77ooo) = ±0.007 as uncertainties related to the -z^ertex 
dependence of the acceptance. 

While the run-to-run variation of the calorimeter energy scale was smaller than 
10~^, the overall energy scale is only known to 10~^ in the 2000 run period ^ . 
By varying the energy scale within ±10~^ we estimated systematic uncertainties 
of ARe(77ooo) — AIm(77ooo) — ±0.001. Furthermore, we investigated the effect of 
possible non-linearities in the shower energy reconstruction by modifying the re- 
constructed energy by an amount AE/E = a/E + I3E + jr, with r being the 
radial distance of the shower from the central detector axis. The allowed ranges 
of the parameters a, (3, and 7 were determined from studies of Kg^, K — > 7r°7r°, 
— > Stt*^, and Tr'^/ry — >■ 77 decays from older data collected in 1999, and found to 
be a = ±10 MeV, /5 = ±2 x 10"^ CeV"\ and 7 = ±10"^ cm^^ When varying the 
constants within these ranges, the fitted value of ?7ooo varied by ARe(?7ooo) = ±0.001 
and AIm(r]ooo) = ±0.002. 

Uncertainties in the acceptance of — > T^^T^^'^'i^m^ decays in the near-target run 
result in an uncertainty of AIm(r7ooo) — ±0.001, while the uncertainty on Re(77ooo) 
is less than 0.001. 

Background processes which could fake a K'^ — > 37r° event play a role only in the 
near-target run, where the majority of decays comes from Ks mesons and A and S° 
hyperons. To estimate the possible background in the Svr^ near-target data 

sample, we compared the tails of the r\„^ distributions of — > 37r° events and 
Ks — > TT^Tr*^ events, which are practically background-free, and found no significant 



In the usual configuration for the direct CP violation measurement, the beginning of 
the decay vohime was defined by a scintillating anti-counter in the near-target beam. In 
2000, this anti-counter was removed to be able to take a 300 times more intense neutral 
beam, and the reconstructed colUmator edge was used as reference position. 
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± 0.001 
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±0.001 


±0.001 


Total: 


±0.015 


±0.017 



Table 1 

Systematic uncertainties on the ryooo measurement. 



discrepancy. Pile-up events could either fake — > 37r° events or lead to accidental 
losses due to the cut on the number of clusters in the calorimeter. Prom studying 
time side-bands and cluster widths, we found no indication for pile-up events. A 
conservative upper limit on the effect of possible background and accidental losses 
in the near-target beam was estimated by loosening the cuts on the energy in the 
hadron calorimeter and the number of clusters in the LKr calorimeter, from which 
we found ARe(77ooo) = ±0.002 and AIm(77ooo) = ±0.009. 

As described in the previous section, the K^K^ dilution was taken from measured 
data of the NA31 experiment. By varying the dilution within the measurement 
errors, we obtained an uncertainty on our fit result of |ARe(?7ooo)| = |AIm(?7ooo)| < 
0.001. 

The systematic uncertainties for r/ooo are summarized in Table 1. The individual 
contributions were added in quadrature, taking correlations between Re(7^ooo) and 
Im(77ooo) into account. 

As a cross-check we have fitted the ratio N'^^^^^^o(t) / N^^^^^o(t) on the far-target 
data only. We obtained Re(77ooo)|far-target = 0.004 ± 0.003 and Im(r7ooo)|far-target = 
—0.003 ± 0.004, with the result being consistent with zero, as expected. 

Two additional independent analyses have been performed, yielding consistent re- 
sults when applying similar cuts. While one of these analyses employed a method 
similar to the one described here [12], the other analysis applied a toy Monte Carlo 
simulation for the geometry correction. 
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Fig. 7. Fit result for ??ooo- The errors include both statistical and systematic uncertainties. 
The lines indicate the exclusion limits for 68%, 90%, and 95% confidence levels. 

5 Results and Discussion 

Our result on r^ooo is 



Re(7^ooo) = -0.002 ± 0.011 (stat.) ± 0.015 (syst. 
Im(ryooo) = -0.003 ± 0.013 (stat.) ± 0.017 (syst. 



(7) 



with a statistical correlation coefficient of 0.77 and an overall correlation coefficient 
of 0.57 between real and imaginary parts (Fig. 7). From this, an upper limit on the 
absolute value of r/ooo is obtained as 



Ir^oool < 0.045 



(8) 



at the 90% confidence level. Using the measured values of Bt{Ki^ Stc^) = 0.211 ± 
0.003 and of the K^^ and lifetimes [2], this result turns into an upper limit on the 
Kg Stt^ branching fraction of 

Br(irs -> 37r°) = \r]ooo\'^ — Br{K^ Sn'^) < 7.4x10"^ (9) 



at a confidence level of 90%. This value is more than one order of magnitude below 
the previous best limit [5]. 

Under the assumptions of CPT invariance Re(77ooo) — Re(e). Fixing Re(?7ooo) = 
Re(e) = 1.66 x lO'^, the fit yields Im(r/ooo)|cpr = 0.000 ± 0.0093t.t ± 0.013,^,t, which 
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[Re(e) - Re(y) - i{Im{x+) + Im((5))] 


-0.001 ± 0.007 


0.005 ± 0.006 


a+_o = 


^ r?;_o Bt{K^ ^ tt+tt-ttO) 


0.000 ± 0.002 


0.000 ± 0.002 


"000 = 


^r?5oo Br(K,^37rO) 


-0.001 ± 0.007 


0.001 ± 0.008 






1.658 lb 0.024 


1.581 lb 0.025 



Table 2 

Parameters «/ = ^(-f^L ^ /)*^(-f^s ^ Df^S for decays of neutral kaons into the final 
state /. r/000 is taken from this measurement. The AS = AQ violating parameter Im(a;+) 
and the CPT violating parameter Re{y) are taken from Ref. [13]. All other measurements 
are from Ref. [2]. 

can be turned into upper limits of |?7ooo|cpr < 0.025 and Bt(Ks — > 37r°)|cpr < 
2.3 X 10"^ at a confidence level of 90%. 

The result on r^ooo fi:"om the fit without imposing CPT can be used to improve the test 
of CPT invariance via the Bell-Steinberger relation [3]. This relation uses unitarity 
to connect the CP violating amplitudes of Kg and K^^ decays with the CP violating 
parameter e and the CPT violating parameter S through 



(l + ^tan0siy)[Re(e)-nm(5)] = Yl A{K, ^ ff A{K, ^ f) /V s 

final 



states / 




final 
states / 



where the super- weak angle is defined by iaxi (psw = 2 (m^ — 7715)/ (F^ — F^). Up to 
now the limit on lm{6) was limited by poor knowledge of r/ooo- Together with available 
measurements of the other CP violating amphtude ratios (see Table 2), the Kj^ and 
Kg lifetimes, branching fractions, and mass difference [2], and the parameters Re(|/) 
and Im(a;+) of possible CPT and AS* = AQ violation [13], the new result presented 
here (Eq. 7) can be used to improve the knowledge of lm[6) and Re(e). Solving the 
complex equation (10) for the free parameters Im{S) and Re(e), taking into account 
correlations between input parameters, yields 

Im{5) = (-0.2 ± 2.0) X 10-^ ^^^^ 
Re(e) = (166.4 ± 1.0) x 10"^ 
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The errors on both Im(5) and Re(e) arc reduced by a factor of 2.5 with respect to 
the previous values [13] and are now hmited by knowledge of 77+_. 

Finally, when assuming CPT invariance in the decay, this can be converted into a 
measurement of the K^K^ mass difference of 

rriK^ - m-^ = (-0.2 ± 2.8) x 10"^^ GeV/c^ (12) 

and yields an upper limit of |m^o— m-^| < 4.7x lO^^^GcV/c^ at the 90% confidence 
level, which is the most precise test of CPT invariance in the mass matrix. 
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